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Abstract 

We quantify the probability per unit Euclidean-time of reversing the magnetization of a vr-Bloch 
vector, which describes the Ferromagnetic Domain Walls of a Ferromagnetic Nanowire at finite- 
temperatures. Our approach, based on Langer's Theory, treats the double sine-Gordon model 
that defines the vr-Bloch vectors via a procedure of nonperturbative renormalization, and uses 
importance sampling methods to minimise the free energy of the system and identify the saddlepoint 
solution corresponding to the reversal probability. We identify that whilst the general solution for 
the free energy minima cannot be expressed in closed form, we can obtain a closed expression for 
the saddlepoint by maximizing the entanglement entropy of the system as a polynomial ring. We 
use this approach to quantify the geometric and non-geometric contributions to the entanglement 
entropy of the Ferromagnetic Nanowire, defined between entangled Ferromagnetic Domain Walls, 
and evaluate the Euclidean-time dependence of the domain wall width and angular momentum 
transfer at the domain walls, which has been recently proposed as a mechanism for Quantum 
Memory Storage. 
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I. INTRODUCTION 



Spin Polarised Scanning Tunnelling Microscopy (SP-STM) measurements of the differen- 
tial conductance profile of Ferromagnetic Nanowires, such as samples of Fe on W(llO) 
highlight the importance of understanding the effective polarization of the surface-probing 
tip subsystem, and its contribution to the measured anisotropy of a Ferromagnetic surface. 
This is because the response of the surface to the magnetic field of the applied probing tip 
is dependent on the local projection of the surface magnetization onto the probing tip, via 
the torque that is induced by the spin-polarized tunneling current. Hence, this mechanism 
is being actively considered for practical Quantum Memory Storage, via the transfer of an- 



gular momentum into Ferromagnetic Domain Walls 



Viemory 



a]. Although the corresponding 



changes in the measured resistivity are not large, it is relatively straight forward to under- 
stand the scattering mechanism induced by thi^ torque in the two limits that the domain 



walls are either very narrow or very broad [3]|8|]. In the first case, spin magnons will be 
specularly reflected without the transfer of angular momentum, whilst in the second case, 
spins adiabatically realign around the local magnetization, whilst precessing with a charac- 
teristic Larmor frequency, and transfer angular momentum into the domain wall. A simple 
(but effective) picture of the local projection of the spins can be parameterized by defining 
a Bloch vector between neighbouring Ferromagnetic Domain (vr-Bloch) Walls, 

sm{ip{x)) = tanh[ (x — Xq) / {w/2) ] , (1) 

where xq is the position of the first domain wall, and w is the domain wall width. This 
relation is obtained by assuming that nonlocal surface effects can be neglected, and by 
minimising the free energy, where the free energy of this system is expressed via the generic 
double sine-Gordon model, 

S = jdx A + MB cos(p + K sinV , (2) 

where A is the exchange interaction stiffness, K the effective anisotropy of the surface, M 
the saturation value of the local magnetization vector, and B the magnitude of the applied 
field which generates Zeeman splitting 9|]. The domain wall width, w, is then obtained 
explicitly by assuming that the dissipative dynamics of the spins obeys the Landau-Lifshitz- 
Gilbert equations, giving. 
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(3) 



K + MB/2 ~ \ V MB 
However, whilst several experiments have confirmed that this locally planar magnon pic- 
ture is essentially correct jll[2|0[6l, this does mean that the domain wall width of the 
Ferromagnetic Nanowire is parameterised in (3) via the non-universal effective anisotropy 
value K. This parameterisation therefore limits the accurate calibration of the angular 
momentum transfer of a Ferromagnetic Nanowire for a practical Quantum Memory Stor- 
age device, and more detailed treatments have stressed the importance of also considering 
noncoUinear magnetization effects [l^ . 

To renormalize this effective anisotropy value, K, we now treat the more general prob- 
lem of minimising nonadiabatic dynamics in (2), and consider the entanglement between 
Ferromagnetic Domain Walls that is induced by rapid changes in the tunneling current. We 
perform a path integral over Euclidean-time for (2) to do this, which is_ made equivalent 



to Imaginary-time in our formalism via an exact Wick-rotation [12] [13|] . Physically, this 
treatment corresponds to the situation where if the system is heated, Stoner excitations 
eventually destroy the ferromagnetic order by overcoming local energy barriers, however, 
rapid recooling is not then guaranteed restore the system to the same initial local valley 
states and so the free energy minima of (2) is not unique [l^. A similar approach to the 



one we present has been developed in 



I ) is not 



to treat the Euclidean-time dependence 



of the nucleation of the Ferromagnetic Domain Walls by applying hanger's Theory. The 
approach is to find the minima of the free energy in (2), parameterized as a dilute gas of 
soliton and anti-soliton states, and to determine the effect of small perturbations to this 
system obtained by adding a single soliton to the grand canonical partition function of the 
system . This gives a probability per unit time of reversing the magnetization of the Bloch 
vector described by the Van't Hoff-Arrehenius law. 
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where Es is the saddlepoint solution for the free energy, (3 is Euclidean-time, n is the 
growth rate of unstable deviations, L is the system length, and {A} are the eigenvalues of 
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the Hamiltonians describing the zero mode, (0), in-plane, and out-of-plane fluctuations 



(p), of the Ferromagnetic Domain Wall Bloch vectors ij]. In order to treat perturba- 



tions of arbitrary size, in our approach, however, we formulate a fugacity expansion for the 
'uU grand canonical partition function, which consists of N soliton and anti-soliton states 
17|[l8|. Minimizing the free energy of this full partition function enables us to normalize 
the probability per unit time in (4) via an additional nonanalytic factor, which corresponds 
to the entropy of adding an additional (A^ — 1) solitons into the system. We then calibrate 
this contribution as a phase factor, via the Wick-rotation of the operators of our formalism 
[ill [3] 13], ^'^is value to quantify the change in the physical scale of the Bloch 

vector in (1) and (3) due to the nonanalytic entanglement contributions (quantum noise) 
arising from magnetization-reversal tunneling transitions between different Ferromagnetic 
Domains across the surface of the Nanowire. 



II. MIXED QUANTUM SPIN CHAIN 



Following [17|, to represent the double sine-Gordon model in (2), we treat the model which 
consists of an antiferromagnetic (AFM) periodic mixed-spin chain of the form, 1 — 1— 3/2 — 
3/2, where 1 and 3/2 represent the spin magnitudes, and the period of the alternation is 
four lattice sites. The Hamiltonian of this model is given by, 

iV/4-l 

TT _ ST T <?(1) -L T <?(^) <?(3/2) J <j(3/2) <-,(3/2) j <j(3/2) ^^(1) 

-H - 2^ -'l,l'^4i •'^4j+l + ^l,3/2»4j+l-^4j+2 + -'3/2 ,3/2'34j+2 •'^4j+3 + -'3/2 ,1 -^4^+3 • -^4^+4 lOJ 

j=0 

where j is the lattice site index of the cell along the spin chain, Ja^b is the nearest neigh- 
bour spin interaction coupling between the neighbouring spins of magnitude a and b, and 
S^""^ is the spin operator of magnitude a. We keep the ratio between like-like and like-dislike 
spin vectors fixed such that Ji^i = ^3/2,3/2 = J and ^1,3/2 = -^3/2,1 = ctJ- The phase space 
of the general quantum spin chain is related to the unstable renormalization group fiow 
of the SU(2) Wess-Zumino-Witten model and the emergence of a gapless groundstate in 
half-integer spin chains - whilst gapped groundstates states are realised by quantum spin 



chains with integer spin values 19] 20|]. These states can be identified in the above Hamil- 
tonian by constructing the symmetric Young tableaux of the spin-3/2 operator (which is 
an SU (2)-spinor) following Haldane's arguments for the destructive interference of quantum 



fluctuations for antisymmetric states 2l| . Tliis analysis reveals that there are three topolog- 
ically distinct realizations of the gapped groundstate of the above mixed-spin model. The 
importance of this particular model is, therefore, that it undergoes two topological transi- 
tions at zero temperature, where the system is driven through the vacuum angle 6* = tt, by 
varying the coupling ratio a. We therefore define the Bloch vectors of this system, not with 
respect to the operator algebra defined between local spin sites, i, but with respect to a path 
integral over Euclidean-time given by the grand canonical partition function, defined by 
summing over all values of the vacuum angle ^, 





r 1 




r rP 






^ V9 exp 


/ drH 
.Jo 


= j exp 


I dr Anicj) 
.Jo 




(7) 



i (n 



Lis 



irs\e) 



{u\e) 



(n|n) 



where and Vn are defined as the compact and noncompact portions of the spin oper- 
ators in (6) determined through the nonperturbative matrix elements A^'^ and (A^'^)', and 
i and r are, respectively, space and Euclidean-time indices [llj[12|. Consequently, there are 
three compact Bloch vectors defined in the regions between the zero temperature transitions 
of the above system, which we use to represent a Ferromagnetic Nanowire consisting of four 
anti-aligned vr-Bloch domain walls. The above lattice operator formalism is explicitly related 
to the effective action of the double sine-Gordon model in (2) via a two-loop expansion of 
the following form. 
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where the Euclidean-time and spatial lattice spacings are quantified via. 
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which quantify the out-of-plane and in-plane fluctuations of the Ferromagnetic Domain 
Walls (Bloch vectors), respectively. 

As indicated by the operator formalism in (9) and (10), since strictly no measurement of 
the inverse spin gap of a finite-size quantum spin chain system can be a divergent function, 
the topological transitions (and hence the Bloch vectors) of the lattice model in (7) and 



(8) are only defined upto the intrinsic lattice IR cutoff scale given by the nonperturbative 
dynamics of A^'^ and (A^"^)'. There is, therefore, a nonanalytic first order region around 
the topological transition value 6' = vr in the quantum spin chain at finite temperatures, 
which represents the mismatch in the overlap of (j) and 6 . Therefore, by finding the 

maxima of (7), we can treat the entanglement of the Ferromagnetic Domain Walls at finite- 
temperatures by quantifying these nonanalytic regions, which correspond to a nonadiabatic 
change in the normalization of </>, from (7) and (8). A similar treatment for renormalizing 
the double sine-Gordon model describing the Ferromagnetic Domain Walls was given in [l^ , 
by treating these nonanalytic regions using a perturbative expansion scheme. 

III. RENORMALIZED ENTANGLEMENT ENTROPY 

We minimize (7), and define an analogous saddlepoint solution to the grand canonical 
partition function treatment in (5), by first evaluating the path integral in (7) using the con- 
ventional continuous-time Quantum Monte Carlo scheme. Practical details of this approach 

n 

are given in [17| , and this procedure is used to generate a set of the nonperturbative matrix 
elements AJ^'^ and (A^"^)' in (8) for the model defined in (6). The important quantities for 
determining the saddlepoint solution in (5), are the eigenvalues, {A(°)}, {A(<^)} and {A^p)}, 
which describe the zero mode, in-plane and out-of-plane fluctuations of the Bloch vectors. 
However, from (10), all of these contributions are incorporated into our nonperturbative 
matrix elements A^*^ and (AJ^'^)'. Therefore, to quantify this saddlepoint solution in our 
new approach we explicitly evaluate the determinant of the operator P, which equivalently 
defines the transfer matrix in the method we use for the importance sampling of the path 
integral. 



where Aj are the eigenvalues of P. The practical reason why only the addition of one extra 
soliton is considered in (5), is that, in general, the support of the eigenvalues in (11) cannot 
be expressed in closed analytic form. However, it is possible to locally regularize (11) via the 
(^-function renormalization prescription |l^. Following this prescription, the logarithm of 
detP in the limit N ^ oo can be defined by analytically continuing the following function, 




(11) 
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from large s to s = 0. Introducing the ^-function, Cp{s) = X]i=i ^r**' logarithm of 
the determinant of P is then completely defined by the first derivative of this (^-function at 
■s = 0. Writing this C-f unction in the usual form of a Mellin transform, 



Cp{s) = YIdI Tr (e-n (13) 

POO 

= / A(^) de (14) 

Jo 

POO 

lndetP= - / A(^) In(^) dO (15) 
Jo 

However, although (15) is a well-defined expression for the path-integral in (7), and the local 
limit of the free energy, because of the rotational symmetry of the Bloch vectors in (8) this 
minima is not unique, and there are several different realizations of the system in (7) which 
are degenerate. For the general solution of (12), not only can the local minima in (15) be 
finite, but also a number of the higher moments corresponding to the tunneling from two-, 
three- and (and more) soliton states. Therefore, to determine the saddlepoint solution of 
(7), we have to evaluate the (^-function prescription at a different point on the fundamental 
strip, at s = 1, for each of these moments. 



/oo 
Aie)Wie)de, j = o...n (16) 

which allows us to address the choices of branch for all of these higher moments that 
relates to the degeneracy of the minima in (15). Collectively minimising these moments, 
at s = 1, then becomes equivalent to analytically continuing the following function via the 
(■-function prescription in (12), 

POO 

H= A{e) \nA{e) de (17) 

Jo 

which corresponds to the familiar Von Neumann form of the entropy of a quantum spin 
system. To find the saddlepoint solution of the full grand canonical partition function in (7), 
it is therefore analytically more well-defined to maximize the entanglement entropy in (17), 
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FIG. 1: The density of the logarithm of the eigenvalues of the transfer matrix in (11), In Aj, plotted 
as a function of the Bloch vector angle (j), defined in (8). 

rather than minimize the free energy in (15), since this resolves the issue of the degeneracy of 
the free-energy minima. However, it follows from the (^-function prescription that the local 
part of the entanglement entropy in (17) is then, in general, a divergent function. Formally, 
the structure of this singularity is a polynomial ring. 

Having maximized the path integral in (7) via importance sampling, we plot the quantity 
InAj as a function of (p, from our definitions in (7) and (11), in Figure 1. This quantity 
defines the normalized density of the logarithm of the eigenvalues along the fundamental 
strip, at s = 1, and corresponds to the entanglement entropy in (17), for the discrete 
system in (7). A more detailed discussion of the system-size of this treatment is presented 
in 17(1, but the general trends evident in Figure 1 (for = 256) are that there is both 
a 0- independent component to InAj, which is relatively fiat and centred on InAj ~ 1.9, 
and two 0-dependent components to InAj, which are asymptotic to ~ and ~ 0.9. 
From (17) these contributions correspond, respectively, to the resolution of the branches 
of the tunneling processes from multiple soliton states at s = 1 (flat plateau), and the 
locally-divergent simple spin- wave magnon component defined in (15) (asymptotes). Via 
the normalization procedure for (p in (8) we can, therefore, use the above numerical value 
of the fiat plateau to quantify the nonanalytic contribution to the Bloch vectors due to 
entanglement between the Ferromagnetic Domain Walls, in the Ferromagnetic Nanowire 
system described by three compact Bloch vectors defined through (6) and (9) by summing 
over all possible tunnelling processes, via importance sampling. 
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IV. ENTANGLED DOMAIN WALL WIDTH 



When the double sine-Gordon model is used to describe the magnons associated with 
the differential conductance profile of ferromagnetic nanowires in jjj^, or more formally 
when this model is used (equivalently) to describe the (infrared) singularities associated 
with topological changes in the ground states of quantum spin chains (where these domain 
walls are viewed as solitons) 19|] 22], in both cases, 0, is used to parameterise the angle 
that describes the orientation of magnetization vectors as they cross the domain walls. The 
two interaction terms in the double sine-Gordon model in (2) therefore have specific generic 
meanings - the first interaction term describes the singlet state of this magnetization vector, 
whilst the second interaction term describes the triplet state found by the Zeeman splitting 
of this magnetization vector. Crucially, when long-ranged dynamical fluctuations are ne- 
glected (in the limits of either very narrow or very broad domain walls), the only physical 
process that can generate this Zeeman splitting is an applied external field. However, when 
these long-ranged dynamical fluctuations are included, this splitting can also be generated 
dynamically by inter-valley tunnelling processes. Thus, to go beyond a simple magnon pic- 
ture for the model requires a treatment of the (nonlocal) dependence of the couplings of the 
two interaction terms, in (2), on the long-ranged dynamical cutoff scale of these inter- valley 
interactions. 

This long-ranged dynamical cutoff scale in our analysis corresponds to relative scale set 
by the flat plateau in Figure 1, and is defined in normalised units of Euclidean time. By 
explicit construction of the operators in (8) this cutoff scale is also defined as a function 
of normalised units of cf). We can infer from our working in the previous section that this 
long-ranged dynamical cutoff scale is at most logarithmically divergent upto the scale of the 
magnon (which also follows directly from the Mermin- Wagner theorem on the absence of 
spontaneous symmetry breaking in < 2 [Ul]). Therefore, in (2), the coupling of the first 
interaction term MB can be smoothly deformed into the couplings of the second interaction 
term K, and the elementary magnons of our analysis our defined by the scale upto which 
they can pass through the domain wall without becoming entangled through inter-valley 
tunnelling. Since in (8) traces out the surface of a sphere in Euclidean time, this smooth 
deformation maps this sphere into an ellipsoid. The eccentricity of this ellipsoid, e, is defined 
from (2) by. 
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FIG. 2: The Lorentz scaling of the effective anisotropy, K, in (1) defined as a function of the 
logarithm of the applied field B using the experimental Ferromagnetic Domain Wall data from Q]. 
The curve indicates the zero entanglement boundary of magnons passing through the domain wall, 
and is obtained from the preceding numerical analysis. 



e = v/l-(r//27 + l)2. (18) 

where '-f = K and rj = MB /2Bmax and the entanglement entropy is defined in units of 
InB = In [B / Bmax) ■ As the elementary magnons of this analysis pass through the domain 
wall they therefore become effectively oblate through entanglement. In the limit that the 
triplet interaction term is larger than the (magnon) singlet interaction term in (2), the 
above eccentricity parameter then becomes purely imaginary - corresponding to a comple tely 
entangled state. A similar elliptical parameterisation of the Bloch vectors is given in 23| . 
to quantify the Pauli cloning of quantum bits. 

In Figure 2 we plot 77/27 + 1 versus logB from the domain wall differential conductance 
profile data obtained experimentally in 2|]. To rescale the effective anisotropy data values 
(which were obtained from fits of the experimental data to the ansatz in (3)) we follow 
the relatively simple prescription of rescaling e, via the graph. We firstly substitute the 
data values of r] and 7 into (18) to calculate the eccentricity, e, of the corresponding Bloch 
vector. We then recalibrate the magnon couplings by defining the new effective anisotropy 
value to be 7' = e, and set the eccentricity of this new Bloch vector to unity by reading 
the corresponding rj' off the graph. Finally, we calculate the rescaled eccentricity value e', 
by putting 7' and rj' in (18). Thus, the ratio e'/e quantifies the entanglement contribution 
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to the effective anisotropy value - as distinct from the Zeeman sphtting. From this (renor- 
mahzation) analysis, the effective anisotropy value differs from the renormalized effective 
anisotropy value by 11% at 200mT, which implies that the level of quantum noise and an- 
gular momentum transfer associated with the domain walls is 11% at 200mT. However, this 
contribution increases dramatically as the field strength is decreased, and magnetic ordering 
is replaced by quantum-mechanical mixing. The fitted line in Figure 2 corresponds to the 
scale-independent numerical value for e we obtained from our numerics in the preceding 
section, which shows good agreement with the data. 



V. SUMMARY 



In this article, we have derived a nonperturbative scaling relation (18) for the angular 
momentum transfer of magnons at ferromagnetic nanowire domain walls from a rigorous 
treatment of the singularities associated with the grand canonical partition function of in- 
teracting Ferromagnetic and Anti-Ferromagnetic Domain Walls. Our scaling relation in (18) 
is defined by finding the infrared cutoff dependence of the couplings in the leading order 



magnon description of changes in resistivity of the domain wall 



17|, and we have used a 



^-function renormalization prescription to parameterise the inter-valley tunneling contribu- 



tions to the magnon couplings via the construction of an exact polynomial ring [13|. The 
advantage of this approach, over existing calculations for the grand canonical partition func- 
tion of interacting Ferromagnetic and Anti- Ferromagnetic Domain Walls [l4|[l^ 16], is that 
the residual degeneracy of the couplings in this polynomial ring means that we do not need 
to mathematically separate the planar and non-planar contributions to the domain wall in 
order to define the magnons, (8). Hence, we are able to quantify the relative changes in 
quantum noise and angular momentum transfer at the domain wall due to the (inter-valley) 
mixing of the magnons. 

Although no approximations are made in this nonperturbative renormalization calcu- 
lation, our analysis rests on the simple magnon picture that is used to parameterise the 
differential conductance profile at the domain wall. We have good evidence in Figure 1 to 
believe that the leading order contribution to the saddlepoint solution of our system is at 
most logarithmically divergent, which is in agreement with the general assumptions of the 
Mermin- Wagner theorem 12] ll[ and the experimental successes of the magnon approach 
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[slj^Q]. However, although we have exphcitly evaluated a grand canonical canonical parti- 
tion function for the interacting Ferromagnetic and Anti-Ferromagnetic Domain Walls con- 
sisting of 256 interacting pairs using importance sampling, the topology of the soliton model 
we have used only contains interaction terms upto third order. Hence, the numerical results 
we have presented may not be consistent with much longer ferromagnetic nanowires, con- 
sisting of a larger number of domains. Nevertheless, our numerical predictions do show good 
agreement with the scaling of the domain wall data in , and our general approach is con- 
sistent with the experimental findings of the leading order magnon behaviour {^[slQljslj?], 
which we are now able to renormalize nonperturbatively. 
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